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ABSTRACT 

The Chandra component of the Great Observatories All-Sky LIRG Survey (GOALS) presently 
contains 44 luminous and ultraluminous infrared galaxies with log (Lir/Lq) — 11.73 — 12.57. Omitting 
15 obvious AGNs, the other galaxies are, on average, underluminous in the 2-10 keV band by 0.7 dex 
at a given far-infrared luminosity, compared to nearby star-forming galaxies with lower star formation 
rates. The integrated spectrum of these hard X-ray quiet galaxies shows strong high-ionization Fe 
K emission (Fexxv at 6.7 keV), which is incompatible with X-ray binaries as its origin. The X-ray 
quietness and the Fe K feature could be explained by hot gas produced in a starburst, provided that the 
accompanying copious emission from high-mass X-ray binaries is somehow suppressed. Alternatively, 
these galaxies may contain deeply embedded supermassive black holes that power the bulk of their 
infrared luminosity and only faint photoionized gas is visible, as seen in some ULIRGs with Compton- 
thick AGN. 

Subject headings: galaxies: active — galaxies: starburst — infrared: galaxies — X-rays: galaxies 



1. INTRODUCTION 

While far-infrared (FIR) luminosity^ is often assumed 
to be a good indicator of the star formation rate (SFR) in 
extragalactic objects (e.g. Young & Scoville 1991, Kenni- 
cutt 1998), it has also been argued that, in the absence of 
active galactic nuclei (AGN), X-ray emission could trace 
the SFR as well; this latter assumption is based primar- 
ily on the observed correlation between X-ray luminosity 
(usually in the 2-10 keV energy range) and FIR luminos- 
ity initially found for nearby actively star-forming galax- 
ies with SFRs -3-30 Mgyr'^ (e.g., RanaUi et al. 2003, 
Grimm et al. 2003; Gilfanov et al. 2004). This rela- 
tion appears to extend to objects with higher SFR (up 
to 100 Mgyr"^) and higher redshift (e.g., Hornschemeier 
et al. 2005; Persic & Raphaeh 2007, Lehmer et al. 2008). 
The hard X-ray emission is generally assumed to have a 
power-law form and to be predominantly due to the col- 
lective emission of high-mass X-ray binaries (HMXBs), 
which are accreting compact objects (neutron stars or 
black holes) formed following the death of short-lived 
massive stars, and therefore expected to have a close re- 
lationship to the starburst. 

There have been, however, indications that objects 
with very high computed SFRs, such as ultraluminous 
infrared galaxies (ULIRGs), might depart from the corre- 
lation because they are underluminous in X-rays for their 
computed SFRs (Persic & Raphaeli 2007, Barger, Cowie 
& Wang 2007). Arp 220, the nearest ULIRG, is a prime 
example, with a 2-lOkeV luminosity 1 dex below the 
general Lx-^fir correlation (see Iwasawa et al. 2005). 
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A further puzzle comes from its X-ray spectrum, which 
shows a strong high-ionization FeK line (mainly Fexxv, 
Iwasawa et al. 2005). This latter finding means that 
Arp 220 is not only X-ray under-luminous, but that the 
hard X-ray emission is not primarily due to X-ray bina- 
ries because of the presence of a high-ionization Fe K line. 
In nearby star-forming galaxies, Fexxv is much weaker 
(e.g., Cappi et al 1999) or undetected, which is consistent 
that the hard X-ray emission is dominated by HMXBs. 
These peculiar properties of the nearest ULIRGs clearly 
warrant further investigation with a larger sample of 
(U)LIRGs. 

The Great Observatory All-sky LIRGs Survey 
(GOALS)^ is a multi-wavelength study of the most lumi- 
nous infrared-selected galaxies in the local Universe, se- 
lected from the flux-limited /i?j4 5 Revised Bright Galaxy 
Sample (RBGS: Sanders et al. 2003). An overview of 
GOALS is given in Armus et al. (2009). C-GOALS (PI: 
D. B. Sanders) is the X-ray component of the project uti- 
lizing data from the Chandra X-ray Observaotry ( Chan- 
dra, hereafter). Details of the X-ray observations are 
described in Iwasawa et al. in prep. Here we report 
results focusing primarily on the Lx- Lyyr, relation and 
the spectral properties of the hard X-ray emission, with 
special attention to the Fe K band. 

The cosmology used to calculate luminosities in this 
paper is = 70 km s'^ Mpc = 0.3, f^A = 0.7, 

based on the latest WMAP results (Hinshaw et al. 2009). 

2. THE SAMPLE 

The current C-GOALS sample is complete down to 
\og{LiY{J Lq) — 11.73, and consists of 44 galaxies from 
the RBGS with redshifts z = 0.010 - 0.088 (see Ta- 
ble 1). We first removed obvious AGN as follows. The 
primary criterion was a flat X-ray spectrum, assessed by 
the X-ray color or hardness ratio (HR). The X-ray color 
is defined as HR = {H - S)/{H + S), where H is the 
2-8 keV counts and S is the 0.5-2 keV counts. Objects 
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TABLE 1 
The C-GOALS sample. 



Object 


logLpiR 


logl-HX 


Object 


logLpiR 


logZ/HX 




Lq 




-1 




Lq 


erg 


HXQ sample 














F17207-0014 


12.42 


41 


34 


NGC 3690 E 


11.49 


41.00 


F19297-0406 


12.39 


41 


26 


ESO 593-IG8 


11.83 


41.28 


P07251-0248 


12.34 


< 40 


90 


VV 705 


11.79 


40.83 


F12112+0305 


12.30 


41 


60 


ESO 255-IG7 


11.75 


41.46 




12 22 


40 


96 




±1. 1 u 


41 15 


F22491-1808 


12.18 


40 


78 


F10173+0828 


11.76 


39.82 


F23365+3604 


12.11 


41 


20 


ESO 203-IGl 


11.84 


< 40.48 


F10565+2448 


11.99 


41 


20 


F01364-1042 


11.82 


41.18 


F15250+3608 


11.98 


< 40 


65 


ESO 239-IG2 


11.73 


40.99 


F09111-1007 


11.97 


41 


11 


P21101+5810 


11.71 


40.30 


ESO 286-IG19 


11.98 


41 


32 


VV 250 


11.68 


41.49 


VII Zw31 


11.89 


41 


43 


F10038-3338 


11.69 


40.80 


ESO 69-IG6 


11.88 


41 


18 


ESO 77-IG14 


11.67 


41.15 


F17132+5313 


11.84 


40 


93 


UGC 4881 


11.65 


40.73 


II Zw96 


11.94 


41 


18 


IC 883 


11.65 


40.81 


AGN sample 














Mrk 231 


12.37 


42 


48 


P19542+1110 


12.06 


42.61 


F14348-1447 


12.35 


41 


77 


ESO 148-IG2 


11.95 


41.92 


P09022-3615 


12.24 


42 


30 


UGC 5101 


11.93 


41.67 


P13120-5453 


12.24 


41 


67 


NGC 3690 W 


11.32 


41.00 


F14378-3651 


12.16 


41 


53 


NGC 6240 


11.81 


42.54 


Mrk 273 


12.15 


42 


40 


ESO 60-IG16 


11.71 


41.86 


F05189-2524 


12.08 


43 


11 


VV 340a 


11.66 


41.46 


F08572+3915 


12.07 


41 


30 









Note. — Source names beginning with " F" or " P" arc from the 
IRAS Faint Source Catalog or Point Source Catalog, respectively. 

with HR > —0.3 are classified as an AGN. This thresh- 
old is chosen because ULIRGs known to host AGN (Mrk 
231, Mrk 273, UGC 5101) cluster just above this value. 
All of the optically identified AGN are selected by this 
criterion. However, Compton-thick AGN are generally 
missed by this criterion because of their weakness in the 
hard band. Therefore, objects that show a strong FeK 
line at 6.4 keV, a characteristic signature of a Compton- 
thick AGN, arc also classified as AGN (NGC 6240, NGC 
3690 West, VV 340a). These criteria classify 15 objects 
as AGN, and they are excluded from further discussion 
leaving an "hard X-ray quiet" (HXQ, as defined by their 
small HR) sample of 30 objects (including NGC 3690 
East). The SFR of these 30 HXQ galaxies, calculated 
assuming that their FIR luminosity is due to dust heated 
by star formation alone, ranges from 60 to 300 Mq yr~^. 

3. L(2-10KEV) AND L(FIR) RELATION 

The sensitivity of Chandra declines steeply above 
7keV. The 2-10 keV luminosity is estimated by extrapo- 
lating the spectral model that describes the data up to 
7keV. When multiple hard X-ray sources are present in 
a single object, those which have no Spitzer-MlPS 24 nm 
counterpart are excluded for calculating the 2-lOkeV lu- 
minosity, assuming they have no relation with the IRAS 
measured luminosity. In NGC 3690, the western and 
eastern galaxies are treated separately due to the differ- 
ence in classification, and their FIR luminosity ratio is 
assumed to be 1:2 based on the 38 /im study (Charman- 
daris et al. 2002). The X-ray luminosity is as observed in 
the rest-frame 2-10 keV band, corrected only for Galactic 
absorption. The median value of the 2-lOkeV luminosity 
is 1.4 X 10^^ erg s~^ (which also matches the mean). 

Figure 1 shows a plot of hard X-ray luminosity, L{2- 
lOkeV), versus FIR luminosity, L(40 — 400/xm), for our 
complete sample of 44 RBGS galaxies listed in Table 1. 
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Fig. 1.— L(2-10kcV) versus L(FIR) for the GOALS AGN (filled 
squares) and HXQ (open circles) galaxies listed in Table 1. For the 
three objects with no 2-10 keV detection, the 95 per cent upper 
limits derived by following Kraft et al. (1991) are shown. The 
dashed line indicates the correlation for the sample of Ranalli et 
al. (2003). The dotted lines show a factor of 2 above and below 
the Ranalli et al. correlation, within which the other studies of the 
same relationship lie. 

For the GOALS HXQ galaxies alone, no clear correlation 

is seen (Kendal's r ~ 0.3) in the limited FIR luminosity 
range. The correlation obtained by Ranalli et al. (2003) 
is re-derived in the same way as done for our sample: log 
L{2 - lOkeV) = -3.84 log L(FIR). This is fully consis- 
tent with the original relation with an offset due to the 
wider FIR wavelength range. The Lx-SFR correlations 
studied by various authors (e.g., Grimm et al. 2003, Per- 
sic & Raphaeli 2007, Lehmer et al. 2008) all lie within a 
factor of ±2 of Ranalli et al. (2003) as indicated in the 
figure. The correlation line is consistent with the upper 
envelope of the HXQ sample. Most of our HXQ galax- 
ies lie well below this correlation. The median value of 
the logarithmic ratio of L(2-10keV) and L(FIR) is -4.5, 
which means our sample of HXQ galaxies is ~0.7dex un- 
derluminous in the 2-lOkeV band for a given L(FIR), or 
corresponding SFR. We note that, given the lack of cor- 
relation in our sample, this is merely a comparison of the 
average value relative to the Lx-.^'Fir correlation. 

4. THE INTEGRATED 4-8 KEV SPECTRUM 

We investigate the integrated hard-band spectrum of 
the 30 HXQ galaxies, since they arc too faint in hard 
X-rays to allow an individual inspection for the iron 
line. Although detection of the FeK line has been re- 
ported for the XMM-Newton spectra of NGC 3690 East 
(Ballo et al. 2004) and Arp 220 (Iwasawa et al. 2005), 
the line docs not have sufficient counts in the Chandra 
data for a significant detection (sec e.g., Clements et al. 
2002). These two objects are included in the integrated 
hard-band spectrum and their contribution to the total 
spectrum is also investigated. In cases where multiple 
hard X-ray sources are detected in a single object (IRAS 
F09111- 1007, ESO 255-IG007, VV 250, ESO 77-IG014, 
IRAS F12112-I-0305, VV 705, UGC 4881), multiple aper- 
tures were selected accordingly. 

Given that our primary interest is in the Fe K prop- 
erties, we restrict tlie analysis to the Fe K and neigh- 
boring band (4.25-7.65 keV in the rest-frame). The FeK 
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Fig. 2.— (a) The integrated 4.25-7.65 kcV spectrum of the 29 
HXQ objects in Table 1. The energy is in the rest-frame. An excess 
in the high-ionization FeK band (primarily due to Fe xxv) at 6.55- 
6.95 keV is clearly visible. The continuum is fiat (F = 1.1 ± 0.5) 
and the FeK band excess has EW ~ 1 keV. (b) The integrated 
spectrum of Arp 220 and NGC 3690 E, two Fe xxv emitters known 
from XMM-Newton observations, c) The integrated spectrum of 
the HXQ sample excluding Arp 220 and NGC 3690 E. Even without 
the two known Fe xxv emitters, a comparable excess flux above 
the continuum in the Fe xxv band is present, indicating that other 
galaxies make a significant contribution to the total line flux in a) . 

complex generally consists of two components with dis- 
tinct ionization states: (a) cold Fe at 6.4 keV; and (b) 
highly ionized FeXXV (6.7 keV) with a minor contribu- 
tion from Fexxvi (7.0keV). These two components can 
be separated well at the spectral resolution of the Chan- 
dra ACIS-S. We define two spectral bands centered on 
the respective line components (6.35-6.55 keV and 6.55- 
6.95 keV), and five neighboring bands for the continuum, 
giving seven spectral bands over the range 4.25-7.65 keV 
in the rest-frame. 

The hard X-ray sources of our sample are generally 
compact, and a small extraction radius, typically ~ 1.5 — 
2", is used. The detected counts have been corrected for 
background using the data from a source- free area on the 
detector in the same observation, although the correction 
is almost negligible with such a small aperture. 

For individual objects, source counts recorded in seven 
bands corresponding to the pre-defined rest-frame bands 
were accumulated. The detector response curve was cor- 
rected by dividing by the mean effective area of each 
spectral band. As most sources have only a few counts 
over the energy range, stacking was done as a straight 
integration of individual sources, neither normalizing by 
the exposure time nor source brightness. 

The total source counts accumulated in the rest-frame 
4.25-7.65 keV band for the HXQ sample are 296 ct. The 
integrated spectrum is shown in Fig. 2a. A strong FeK 
line is immediately recognized with a 3a excess. The line 
peaks at the high-ionization FeK band (6.55-6.95 keV), 
indicating that the line is primarily due to Fe xxv from a 
highly ionized medium. The total counts collected in this 
line band is 35 ct. There might also be a slight excess in 
the 6.4 keV line band, which is however not significant. 

Given the limited statistics, a further analysis beyond 
the line detection is not warranted. However, given our 
stacking method, which is biased for sources that are 
bright and/or with a long exposure, it is at least worth 



examining whether Arp 220 and NGC 3690 E, the two 
known Fe K emitters from their XMM-Newton spectra, 
dominate the line detection. These two objects and the 
rest of the HXQ sample were integrated separately, and 
their spectra are presented in Fig. 2b and 2c. Both spec- 
tra show a ~ 2(T excess in the Fe xxv band with compa- 
rable line fluxes. This means that, besides Arp 220 and 
NGC 3690 E, there is a significant contribution by other 
galaxies to the detected line flux. The higher continuum 
level in the spectrum of the sample without Arp 220 and 
NGC 3690 E (Fig. 2c) suggests that some sources without 
strong FeK emission are present, and that they may be 
dominated by HMXBs as found in nearby star-forming 
galaxies. 

Fitting a power-law to the continuum gives a photon 
index of F = 1.1±0.5, indicating a rather hard conitnuum 

above 4keV. The equivalent width of the excess in the 
Fe K band with respect to the power law continuum is 
EW = 0.9 ± 0.3 keV. The 90% upper limit of the EW 
for a 6.4 keV line is 120 eV. 

Fitting a thermal spectrum (MEKAL) gives a temper- 
ature of kT = 8 ± 2 keV and an Fe metallicity consistent 
with the solar value. Considering that the integrated 
spectrum is diluted by some sources that contribute only 
to the continuum emission (see above), the Fe metallicity 
of the line emitting sources is likely to be super-solar. 

Regardless of the choice of spectral model, the strong 
high-ionization Fe K line is not compatible with the X- 
ray spectrum of known HMXBs. FeK emission is of- 
ten observed from HMXBs in our Galaxy, but the major 
component is the cold line at 6.4 keV and its EW is on 
average -200-300 eV (e.g. White et al. 1983). The much 
stronger, high-ionization line in the integrated spectrum 
of the HXQ sample suggests that X-ray binaries are not 
the primary source of the 4-8 keV emission. 

5. DISCUSSION 

The general quietness of hard X-ray emission in the 
HXQ galaxies (Fig. 1) is likely related to the detection of 
the high-ionization Fe K feature, which rules out HMXBs 
as the primary origin of the faint hard X-ray emission. 

The powerful FIR emission in the HXQ (U)LIRGs pre- 
dicts an abundance of HMXBs if star formation is the 
dominat power source. The Chandra spectra, however, 
suggests these sources contribute little to the hard X-ray 
emission. Therefore, either the bulk of the HMXBs are 
missing or they are heavily obscured from view. Suppose 
heavy obscuration is the cause. Since the Fe K band is 
still transparent to A'h'^ 10^^ cm~^, the obscuring col- 
umn has to be larger, perhaps the order of 10^^ cm~^. 
A dense, nuclear molecular disk (with a dynamical mass 
of the order of lO^M©) is often found in the central part 
of ULIRGs (e.g., Bryant & Scoville 1999) and thus, such 
a large column may not be unrealistic. The size of the 
nuclear molecular disk is generally found to be a few 100 
pc but the densest part can be as small as 30 pc, as found 
by high-resolution studies of Arp 220 (Downes & Eckart 
2007, Sakamoto et al. 2008, Aaho et al. 2009), where 
a rough estimate of the density is lO'^-lO® cm^''. Ac- 
cording to Grimm et al. (2003), the number of HMXBs 
with Lx > 2 X 10^^ ergs~^, which would dominate the 
integrated luminosity of the entire HMXB population, is 
~ 3OO(5Fi?/lOOM0 yr-i). For SFR = 200Mq yr^^, as 
estimated for Arp 220, the number is — 600, and most 
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of them need to be confined within 30 pc in order for 
their radiation to be suppressed. The required stellar 
density is unusually high and could be a problem in this 
particular case. The other galaxies may have a larger 
starburst region and could avoid this problem provided 
the required mass of the obscuring gas does not exceed 
the dynamical mass. 

Once the HMXBs are obscured from view, the Fe xxv 
spectrum can be explained by high-temperature (T ~ 
10* K), thermal gas produced by a starburst. Possible 
sources are (i) an internally shocked hot bubble* pro- 
duced by thermalizing the energy of supernovae (SNe) 
and stellar winds, as predicted by, e.g., Chevaher & Clegg 
(1985); and (ii) collective luminous SNe. The latter is 
ruled out if HMXBs are all embedded in heavy obscu- 
ration because the SNe should also be embedded in the 
same obscuring material. The former would excavate 
the obscuration and become visible. With the assumed 
high SFR, the luminosity and the spectrum with strong 
Fexxv can be reproduced (e.g., Iwasawa et al. 2005 
for Arp 220). The diffuse hard X-ray emission seen in 
M 82 (Griffiths et al. 2000, Strickland & Heckman 2007, 
Ranalli et al. 2008), after removing resolved X-ray bina- 
ries, has, in fact, a comparable X-ray to FIR ratio with 
that of our GOALS HXQ sample. The strong Fe line is 
also well matched with the prediction for metal enriched 
gas of this origin. 

With the presence of dense nuclear g heavily ob- 
scured AGN is also a possible explanation for the X- 
ray spectra of our HXQ sample. Massive black holes in 
the process of rapid growth are naturally expected to be 
present in (U)LIRGs. When a central AGN is deeply 
buried with a covering factor nearly unity, reprocessed 
light from the obscuring matter would have difficulty es- 
caping. The high-ionization FeK line could then origi- 
nate from extended, low density gas which is photoion- 
ized by the AGN. In a number of Compton thick Seyfert 

^ This should be distinguished from a much more extended, soft 
X-ray nebula, which most likely originate from the swept-up inter- 
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2 galaxies, e.g., NGC 1068, Fexxv has been seen (but 
normally weaker than the cold line at 6.4 keV). This is 
in fact the feature originally predicted by Krolik & Kali- 
man (1987) when assuming a likely scattering medium 
to produce the polarized broad-line region in NGC 1068 
(Antonucci & Miller 1985). Where the cold X-ray re- 
flection is suppressed, the weak, highly photoionised gas 
could be the only observable feature. We note that 
Fexxv has been observed as the primary FeK feature 
in some Compton-thick AGN residing in the ULIRGs 
IRAS F00183-7113 (Nandra & Iwasawa 2007, Ruiz et 
al. 2007), The Superantennae (Braito et al. 2009), and 
possibly UGC 5101. 

In Arp 220, recent measurements of a compact dust 
emission source with a steep temperature gradient in 
the western nucleus provides an argument for an AGN 
(Downes & Eckart 2007, Aalto et al. 2009; but see 
Sakamoto et al. 2008). Millimeter wavelength molecular 
line observations also suggest the presence of X-ray dom- 
inated chemistry, which favors heating by AGN (Aalto et 
al. 2007, Imanishi et al. 2007; e.g., Meijerink & Spaans 
2005 for theory). The similarity of the X-ray properties 
suggests that the same explanation might apply for at 
least some of the GOALS HXQ galaxies - i.e., despite 
the lack of outward evidence of an AGN, a significant 
fraction of the infrared output could be powered by a 
heavily obscured AGN. 
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